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Abstract. We present ion beam erosion experiments performed in ultra high vacuum 
using a differentially pumped ion source and taking care that the ion beam hits the 
Si(OOl) sample only. Under these conditions no ion beam patterns form on Si for angles 
d < 45° with respect to the global surface normal using 2 keV Kr + and fluences of 
ss 2 x 10 22 ions/m 2 . In fact, the ion beam induces a smoothening of preformed patterns. 
Simultaneous sputter deposition of stainless steel in this angular range creates a 
variety of patterns, similar to those previously ascribed to clean ion beam induced 
destabilization of the surface profile. Only for grazing incidence with 60° < i? < 83° 
pronounced ion beam patterns form. It appears that the angular dependent stability 
of Si(OOl) against pattern formation under clean ion beam erosion conditions is related 
to the angular dependence of the sputtering yield, and not primarily to a curvature 
dependent yield as invoked frequently in continuum theory models. 
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1. Introduction 

Ion beam surface patterning of amorphous and crystalline materials has attracted 
considerable interest in the recent years. This interest is of twofold origin. First it 
results from the intricate physics involved in ion beam patterning which has not yet 
been entirely uncovered. Second, ion beam nanopatterning is a relatively cheap and 
low tech method with a number of potential applications ranging from an anti-reflection 
surface finish (TJ El [3] over orienting large molecules [I] to nanomagnetism [5_1 EJ [7J E] . 
As Si is a prime material of technology and readily available in high purity and quality, 
it is not surprising that the ion beam patterning studies are numerous for Si. At room 
temperature Si readily amorphizes during ion exposure [9]. The loss of anisotropy 
and crystal structure appeared to make it an ideal material to be described by the 
continuum theory of ion erosion, which effectively averages out atomistic details of the 
processes. Early theoretical work in the continuum theory approximation considered the 
dependence of the sputtering yield Y on the angle 9 of the ion beam with the respect to 
the local surface normal as a decisive factor for surface morphological evolution ([rO|[TT] 
and references therein). Note that we distinguish here and in the following between the 
local angle of incidence 9 (measured with respect to the surface normal of a specific 
surface element) and the global angle of incidence *& (measured with respect to the 
normal of the average surface plane). Since the seminal publication of Bradley and 
Harper [12] the dependence of the sputtering on surface curvature was considered to 
be the key for pattern formation. While there are certain conditions, where ion beam 
erosion of Si does not cause pattern formation [131 EL the overwhelming number of 
investigations find pattern formation on Si in a large parameter space. One class of 
prototypical patterns are dot or hole patterns observed for normal incidence noble gas 
ion erosion with energies up to a few keV and at temperatures in the amorphization 
regime [HI [151 El El El E] ■ Another class of exemplary patterns are ripple ones with 
the ripple wave vector k parallel to the ion beam azimuth. Quite some observations here 
refer to elevated temperatures in the crystalline regime [2"0"l |2"T1 12"2"] . However, also at 
room temperature in the amorphous regime ripple patterns were observed [231 EU ESI EE] , 
the most regular ones for slight off-normal conditions with d f» 15° using noble gas ions 
with energies up to a few keV [HI [271 EH1 EH]. The variety and complexity of observed 
patterns stimulated the development of continuum theories [301 EH EH EH El] extending 
the concept of the curvature dependent yield as prime destabilization mechanism. 

Co-deposition of trace amounts of foreign species during ion erosion has been found 
already long time ago to give rise to microstructure formation [321 ESI EH EH]- It is 
currently used as a tool for surface texturing [391 HO]- Recently it became obvious 
that impurities may influence pattern formation unintentionally. Mo co-sputtered from 
sample clips during ion erosion was found to foster dot formation at normal incidence 
[HI H21 H3]. The ion flux and fluence were found to affect the amount of Mo and Fe 
deposited on the sample from the ion source and thereby to change normal incidence 
patterns [lU H5] . 
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In this situation the question arises, whether there is a hidden chemical or impurity 
factor in pattern formation on Si. Such a largely disregarded attribute could make it 
impossible for theory to come up with an adequate material parameter based description 
of pattern formation. We therefore conducted an erosion study under clean surface 
science conditions, using ion beam parameters for which we expected pattern formation 
to take place. 

2. Experimental 

The experiments were performed in a variable temperature scanning tunneling 
microscopy (STM) apparatus [IB] with a base pressure < 6 x 1CT 11 mbar equipped with 
a differentially pumped ion source, a low energy electron diffraction (LEED) system, and 
a load-lock for Si sample transfer. For erosion the samples were exposed to a 2 keV Kr + 
fine focus ion beam with full width at half maximum « 1 mm scanned over a sample area 
of about 4 mm x 4 mm resulting in a time averaged ion flux of 5 x 10 16 ionsm~ 2 s _1 at 
300 K at angles 0° < d < 83°. The fluence was F« 2x 10 22 ionsm" 2 . Due to differential 
pumping during ion exposure the working pressure was below 9 x 10 -8 mbar. The ion 
current onto the sample was controlled with a Faraday cup, which could be moved 
precisely into the sample position. After ion exposure the pressure dropped quickly 
into the 10 -11 mbar range and imaging by STM was performed subsequently. For the 
co-sputter deposition experiments a piece of target material was mounted vertically on 
the sample. The ion beam was impinging at d = 30° onto the sample surface and with 
an angle "d^ep = 60° onto the sputter target. In the co-sputter deposition experiments 
the fluence was only F ~ 5 x 10 21 ionsm~ 2 , but still sufficient to guarantee pattern 
development. Again the samples were analyzed in situ by STM and subsequently 
analyzed ex situ by atomic force microscopy (AFM), secondary electron microscopy 
(SEM), a compact phase-shifting interferometer and secondary ion mass spectrometry 
(SIMS). For SIMS we used 250 eV for depth profiling. Quantitative image analysis 
was conducted by WSxM [47] and differential sputter yield and energy distribution 
calculations were performed with TRIM.SP |48j . 

3. Angle dependent smoothening and pattern formation on Si(OOl) 

The -^-dependence of the morphology after ion bombardment is shown in Fig. [T] 
Unexpectedly, for 0° < $ < 45° no patterns form and the root mean square roughness 
a remains very low, a ~ 0.2 nm [compare Figs. [l]^a) - pud)]. For $ = 60° as shown in 
Fig. [lje) ripples with a small amplitude and wave vector k parallel to the ion beam 
azimuth developed. From the power spectral density of large topographs we obtain an 
average wavelength of A ~ 46 nm. Still the surface does not destabilize to any significant 
extent: a is just 0.5 nm after removal of the order of 1 /im of material. Pronounced 
pattern formation takes place in a narrow angular range 75° < °d < 79°, just around the 
angular range of maximum sputter yield Y (■#). At $ = 75° represented by Fig. [IVf) we 
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Figure 1. STM topographs of Si(100) after a fluence of 2 w 10 22 ionsm- 2 2keV Kr+ 
at 300 K. The angle of incidence I? with respect to the surface normal is (a) 0°, (b) 
15°, (c) 26°, (d) 45°, (e) 60°, (f) 75°, (g) 79° and (h) 83°. The image size is for (a)-(d) 
and (h) 316 nm x 316 nm and in (e)-(g) 625 nm x 625 nm. The white arrows in (b)-(h) 
indicate the ion beam azimuth. The corrugation Az is 3nm in (a)-(e) and (h) and 
30 nm in (f) and (g). Inset in (f): height profile along line indicated in (f). 



find a pronounced sawtooth profile ripple pattern with a = 6.4 nm. The sawtooth profile 
displays extended facets forming a small angle of ~ 7° with respect to the global surface 
plane resulting in an angle 9 m 82° of the ion beam with respect to the local surface 
normal [compare inset of[IJf)]. The smaller facets form an apparent angle of ~ 19° with 
respect to the surface plane, resulting in an apparent angle 9 w 56° of the ion beam 
with respect to the local surface normal. Due to convolution of the surface profile with 
the STM tip of finite sharpness the apparent 9 ~ 56° is just an upper bound to the true 
angle of the ion beam with respect to the local surface normal, which is likely to be much 
lower. At $ — 79° represented by Fig. [T](g) instead of ripples we find roof tile structures 
[30] ■ The structures are now elongated along the ion beam and if one would like to 
assign a A; to them (which is not justified), it would now be normal rather than parallel 
to the ion beam azimuth. The roughness is with a = 5.3 nm similar to the d = 75° 
case. Note also that the facet structure of the roof tiles in Fig. [T](g) has similarities, 
but also discrepancies compared to Fig. [ljf). Also the roof tiles display extended facets 
parallel to the ion beam direction with a small angle ~ 3° with respect to the global 
surface plane resulting in 9 pa 82°. The facets normal to the ion beam direction visible 
in Fig. [l](f) became arrow tips. The angle of the arrow tips with respect to the average 
surface plane is ~ 19° limited again by the surface profile - tip convolution. In Fig. [l](h) 
for $ = 83° the surface is extremely smooth again with a = 0.2 nm. A faint ripple 
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Figure 2. (a) Roughness a as a function of the angle of incidence $ for the experiments 
represented by Fig. [l] (b) Sputtering yield Y(9) as a function of the angle 9 of the ion 
beam with respect to the local surface normal for 2 keV Kr + on Si as calculated by 
TRIM.SP @g. Lines to guide the eye. 

— * 

pattern with k normal to the ion beam azimuth is visible; however, the amplitude of 
the pattern is marginal. If forced to define a critical angle of ripple rotation we would 
set it to d c = 77°, in reasonable agreement with previous work [40J finding a $ c = 80° 
for 5 keV Xe + ion erosion of Si. 

Figure [2]^a) summarizes the evolution of a already discussed during the presentation 
of the STM data above. The absence of roughening for d < 45° is no accident, but due 
to an inherent smoothening action of the ion beam. To support this idea after STM 
imaging of the rough d = 75° ripple pattern with a = 6.4 nm, we exposed it to a fluence 
of ~ 5 x 10 21 ionsm -2 of 4keV Kr + at normal incidence ($ = 0°) and halved thereby 
the roughness. 

In order to discuss the mechanism of pattern formation for 60° < $ < 79° we 
consider the angular dependence of the sputter yield Y on the local angle of incidence 
9 as calculated by TRIM.SP and shown in Fig. [2^b). Evidently, this distinction is only 
relevant, when the surface is not flat. Y(9) displays a broad local minimum for 9 0°, 
an absolute minimum for 9 = 90° with Y(90°) = and a strong maximum around 
9p w 75° (see also [19]). It is obvious, that pronounced patterns develop only for global 
angles of incidence d close to the maximum of Y, i.e. where the erosion rate is largest 
and where the erosion rate depends strongly on This points to the application of 
the deterministic approach of erosion profile evolution, where the sputtering yield is 
assumed to be a function only of the local surface gradient [SB ED QUI ITT] . 

If the sputter yield depends only on the local angle of incidence 9, a smooth surface 
eroded under a global angle can become unstable, if fluctuations are present which give 
rise to locally varying sputtering yields. An analysis of surface topography evolution 
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assuming that the sputter yield only depends on the local incidence angle 9 shows that, 
due to fluctuations, an initially flat surface can decompose into a faceted profile which 
displays facets where the sputter yield is a minimum or maximum [501 Elj - For these 
facets which predominate during ion erosion the local ion incidence angles 9 are given 
by 0°, ±0p, or 90°, respectively. It should be noted, that 90° facets can only form if in 
the initial surface profile ion incidence angles 9 > 9p are already existing (501 E2]- For 
surfaces with small height fluctuation this condition is only achieved if •& is close to 75°. 

As the average surface orientation is maintained during erosion, in the angular range 
d pa 75° — 80° the surface profile decomposes into patches with local ion impingement 
angles 9 close to the minima of Y(6), i. e. into facets 1 and 2 with 9\ < d < 2 - 
Considering the example represented by Fig. [l](f) with •§ = 75° the flat surface 
decomposes into a sawtooth profile with facets forming a large angle with the ion beam 
9 1 < 54° and facets nearly parallel to the ion beam (92 ~ 82°). For both 9\ and 92 the 
yield Y is much smaller than for 9 = 75° [compare Fig. j2^b)]. It is evident that such 
a sawtooth profile considerably reduces the global erosion rate. However, the observed 
local angles 9 = 54° and 9 = 82° do not match the theoretical predictions of 9 = 0° 
and 9 = 90°. For the steep facet we attribute this difference largely to inability of our 
STM tip to measure the proper angle. For the facet with the small slope the difference 
of the measured 82° and predicted 90° can not rely on a measurement problem. The 
measured difference may have the following reasons: (i) The fluence used was too low to 
allow the surface to reach the dynamic equilibrium and the ensuing facets; (ii) There is a 
significant uncertainty in the TRIM.SP calculations of Y(9) for very grazing angles due 
to the fast change of the yield with 9 and its sensitive dependence on surface structure. 
The yield at 9 = 82° might thus be already close to zero, such that the driving force 
(minimization of erosion rate) has largely vanished and a further change of the facet 
angle is kinetically frozen. 

In contrast, for near normal ion incidence angles d ~ 0°, 9 < 9p and surfaces 
with small height fluctuations are stable because profiles with 9 = 0° are steady state 
surface configurations [SHI EE]- Furthermore, as hillocks are eroded faster than valleys 
(contrary to the Bradley-Harper theory [12]) the ion beam induces a smoothening of 
rough surfaces in this angular range as observed by us experimentally (compare also 
Fig. 4 of p] and [53]). 

For the case of room temperature ion beam erosion of Si with 2keV Kr + , the 
angular dependence of Y(9) appears to explain the ranges of stability observed as well 
as the orientation of the pattern facets formed in the unstable regime. So far we left 
open, which kind of fluctuations might initiate the faceting of the surface for $ pa 9 p. 
While the stochasticity of the ion impacts gives rise to fluctuations by itself, it may 
be that other effects contribute to fluctuations, e.g. a curvature dependence of the 
sputtering yield. However, the pattern formation scenario we observe here can hardly 
be reconciled with models based on Bradley-Harper theory. A curvature dependent yield 
as destabilization mechanism neither explains the extended angular range of stability nor 
the rather abrupt angle-dependent transitions from smooth surfaces to faceted patterns 
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and back to smooth surfaces, i.e. abrupt transitions from surfaces with zero curvature 
to ones with curvature singularities. By invoking the angular dependence of Y(6) to 
explain surface destabilization, no predictions related to time scales necessary for pattern 
evolution nor to characteristic length scales in the pattern are possible. We therefore 
consider the angular dependence of Y(9) only as one important element in a theory of 
pattern formation on Si. Secondary effects - e.g. surface diffusion, ballistic drift, viscous 
flow, etc. - must probably be considered to be responsible for selection of characteristic 
scales. 

With ion energies of the order of IkeV and for noble gases like Ar + , Kr + or 
Xe + and angles i? < 30° a large variety of dot and ripple patterns was obtained by 
us pH EH EHl [29] and other groups pH EH HZl US 122] ■ We consider these results 
to be at variance with the present findings. It was noticed recently that the relative 
concentrations of Mo and Fe emerging from the ion source during sputtering influence 
the pattern appearance at $ = 0° [HJ US]. However, whether patterns evolve at 
all without impurities was not considered. Our observations are consistent with the 
observation of Mo-seeding at $ = 0° [4T"j W2\ H3] and the absence of patterns without Mo- 
seeding [HI H2]. We claim here that patterns are entirely absent without intentional or 
unintentional co-deposition of impurities after noble gas bombardment on Si(001) with 
energies of the order of IkeV and for -d < 45°. Below we provide additional evidence 
for this statement. 

4. Angle dependent pattern formation on Si(OOl) with co-sputter 
deposition 

To substantiate our claim that impurities resulting from the sputtering process are in 
fact responsible for the great number of patterns resulting for ion erosion with $ < 45° 
we performed dedicated co-sputter deposition experiments. As sketched in Fig. |3](a) the 
ion beam was impinging at $ = 30° onto the sample surface, i. e. at an angle where 
no pattern formation is expected. Additionally the ion beam hits a piece of stainless 
steel (Fe 84% and Cr 13%) at an angle $dep — 60°, from which some material is sputter 
deposited onto the eroding Si surface. This setup is similar to the arrangement used 
for surfactant sputtering [10] . The STM topographs were taken along a line normal to 
the center of the steel plate. As obvious from Figs. [3]^b)-(g) the resulting morphologies 
strongly depend on the normal distance x to the stainless steel plate. Fig. ^b) displays a 
relatively smooth surface with small hole structures as highlighted by the inset (compare 
also [44J ) . With increasing x the roughness a shoots up beyond 10 nm and ripples with 
k parallel to the ion beam azimuth form [compare Figs, [^c) and (d)]. Upon further 
increase of x the roughness a decreases again and the ripple pattern transforms [Fig.pTe)] 



to a dot pattern [Fig. pff)]. Eventually, for the largest distance measured a is comparable 




to a situation without co-sputter deposition and patterns are absent [compare insets of 



Not every co-sputtered material induces patterns. We used the same set-up with a 




Figs.gb) and (g)]. 
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Figure 3. (a) Schematic sketch of co-sputter deposition geometry. The recording 
positions of the STM topographs shown in (b)-(g) are indicated, (b)-(g) STM 
topographs of Si(OOl) after a fluence of « 5 x 10 21 ionsm~ 2 2keV Kr+ at 300 K at 
d = 30°. The image width is always 625 nm; the inset size in (b) and (g) is 100 nm x 
100 nm The ion beam azimuth is indicated by a white arrow. The corrugation Az is 
(b) 2 nm, (c) 45 nm, (d) 20 nm, (e) 8nm, (f) 4nm and (g) 2nm. 




x (mm) 

Figure 4. Roughness a (open dots, left y-axis) and feature separation A (open 
triangles, right inner y-axis) of morphologies obtained in the co-sputter deposition 
visualized in Fig. [3] as a function of the distance x from the stainless steel plate. The 
Fe concentration cp c normalized to the Si concentration as function of x is shown on 
the outer right y-axis (open squares). Lines to guide the eye. 

plate of highly oriented pyrolytic graphite instead of stainless steel. The entire Si wafer 
remained smooth with a roughness o « 0.2 nm and no patterns formed. 

Figure [4] displays the roughness a (open dots, left y-axis), the feature separation 
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Figure 5. (a) Top view sketch indicating the positions, where images (b)-(d) were 
taken by scanning electron microscopy (SEM) after the erosion experiment visualized 
by Fig. [3] had been conducted, ir-axis indicates the normal to the center of the steel 
plate. The direction of e-beam illumination is at 45° with respect to the surface normal 
and the e-beam energy 15keV. The image width in (b)-(d) is 2 fim. 

A (open triangles, right inner y-axis) and the Fe concentration cp e (open squares, right 
outer y-axis) for the patterns on the Si wafer as function of x. Schematically, four 
pattern ranges are distinguished in Fig. |4} nanoholes with x = Omm — 1.0 mm, ripples 
with x = 1.0 mm — 2.5 mm, dots with x = 2.5 mm — 3.8 mm and a smooth surface 
for x > 3.8 mm. The roughness displays a pronounced maximum in the ripple range 
with a w 10 nm, fades gradually away in the dot range with a < 2nm and is very 
low in the nanohole and smooth surface range with o ~ 0.3 nm. The characteristic 
feature separation varies unspectacularly from a typical ripple wavelength A f» 80 nm 
to a typical dot separation with A m 50 nm. Based on the variety of patterns we 
expected an inhomogeneous Fe concentration cp e . As Rutherford backscattering needs 
a much larger sample spot (a few mm) for chemical analysis, we decided to perform 
secondary ion mass spectrometry. Quantification of our data is based on the tables of 
SIMS impurity signals in Si obtained for 8keV sputtering (5U [55] . Here we used 
250 eV to mimimize implantation and mixing and achieve a high depth resolution. 
The surface roughness and the potentially inhomogeneous Fe distribution makes our 
quantification problematic, which therefore should be considered with caution. We 
measure a maximum of the relative iron concentration cp e ~ 0.13 in the ripple range. 
Still cfc ~ 0.02 is sufficient to induce the dot range. It appears plausible that due to ion 
beam induced mixing the co-sputtered Fe forms a silicide in the surface near layers as 
it has been found for normal incidence ion bombardment through an alternating cold 
cathode ion source [44j. 

We further explored the pattern at locations displaced from the normal to the 
center of the stainless steel plate (i.e. left and right of the x-axis) by scanning electron 
microscopy. Figure [5ta) sketches the locations of the SEM micrographs displayed in 
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Figjijb)-(d). It is obvious that the ripple wave vector k depends on the position with 
respect to the stainless steel plate. Apparently k is parallel to the average azimuthal 
direction of the atoms sputter deposited from the stainless steel plate onto the Si wafer. 

The co-sputter situation is complex and in order to prepare an understanding of 
pattern formation in this situation it is mandatory to consider energy, direction and 
species of the particles impinging onto the Si wafer. A given location on the Si wafer in 
front of the stainless steel plate receives 2 keV Kr + ions, sputtered Fe atoms (we neglect 
Cr atoms in the following discussion for simplicity) and Kr particles reflected from 
the stainless plate. Fig. J6^a) displays the differential Fe sputter yield ^fK^) averaged 
over the azimuthal angular range tp G [—30°, 30°] as calculated by TRIM.SP [38]. The 
distribution of the Fe atoms is anisotropic with a broad peak of emission centered at 
about $ pa 37°. As the Si wafer can be imagined horizontally under the polar diagram, 
the main flux of sputtered Fe atoms is directed towards the Si wafer. The total yield 
Ype P er 2 keV Kr + is substantial and amounts to Y Fe = 6.4. 

Fig. 6^b) shows the differential backscattering coefficient ^f 31 ^) of the Kr + 
particles backscattered from the stainless steel plate and averaged over the azimuthal 
angular range (p G [—30°, 30°]. The distribution of the backscattered particles is peaked 
at an exit angle of $ pa 68°. In average the backscattered Kr particles will hit the Si 
wafer closer to the stainless steel plate than the sputtered Fe atoms, i.e. at a smaller 
distance x. The total backscattering coefficient amounts to Rkt = 0.24. 

Fig. [6^c) displays the energies E of the sputtered Fe atoms and of the backscattered 
Kr particles as a function of the polar emission angle ■&. The data plotted is averaged 
over the azimuthal angular range p> G [—30°, 30°]. The energy of the sputtered Fe atoms 
is moderate with a maximum of ~ 210 eV for d pa 74°. Note that at the most probable 
angle of emission of sputtered Fe atoms $ ~ 37° the energy of the sputtered Fe atoms is 
only about 65 eV in the presented azimuthal angular range. The average energy of the 
Fe atoms sputtered in the azimuthal angular range ip G [—30°, 30°] is 112 eV, the average 
energy of all atoms sputtered in forward direction with ip G [—90°, 90°] is 72 eV. Though 
substantial, the energy of the Fe particles is insufficient to cause significant sputtering 
of the Si wafer. The energy of the backscattered Kr particles increases monotonically 
with their emission angle *& and reaches ~ 1200 eV for •§ pa 90°. Note that at the most 
probable emission angle of i? pa 68° the energy of the emitted particles is about 755 eV. 
The backscattered particles will therefore cause strong erosion of the Si wafer near 
the steel plate. The average energy of the backscattered Kr particles in the presented 
azimuthal angular range ip G [—30°, 30°] is 583 eV, the average energy of all particles 
backscattered towards the Si plate in the angular range (p G [—90°, 90°] is 470 eV. 

Fig. |7[a) displays a section of a height topograph of the Si-wafer after co-sputter 
deposition measured by phase-shifting interferometry. We attribute the hill in the height 
profile in the range 1mm < x < 2.5 mm with a height of ~ 75 nm to the significant 
amount of co-sputtered Fe in this range. Note that the hill does not consist entirely of 
Fe - actually most of the co-sputter deposited Fe has been resputtered by the Kr + ion 
beam. However, the co-sputter deposited Fe partly shielded the Si from being removed 
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Figure 6. (a) Differential sputtering yield ^jff- and (b) differential backscattering 
coefficient ^ of 2 keV Kr+ impinging on the steel plate plotted as a function of 
the polar emission angle $. The data shown is averaged over the azimuthal angular 
range p e [—30°, 30°]. The arrows indicate the direction of the impinging Kr + . (c) 
Energy of the sputtered Fe atoms or backscattered Kr particles as a function of their 
polar emission angle The data shown is averaged over the azimuthal angular range 
p e [-30°, 30°]. All data calculated by TRIM.SP gS]. 



by the direct Kr + ion beam resulting in a significantly reduced erosion depth. The 
trough next to the stainless steel plate for < x < 1mm results from the additional 
sputtering caused by the energetic component of the Kr particles reflected from the 
stainless steel plate. According to mirco phase interference patterns at the edge of 
the scanning area of the ion beam the average erosion depth is ~ 300 nm. The SEM 
topograph in Fig. [7](b) shows a contrast change at the onset of the co-sputter hill caused 
by the transition from the smooth nanohole to the rippled surface. This is also consistent 
with the pattern sequence as observed by STM in Fig. |4j In view of the significant flux 
and energy of the reflected Kr particles arriving in the trough area, there the lower cp e is 
a consequence of the enhanced resputtering. Dot patterns form beyond the co-sputter 
deposition hill, where the arriving flux of Fe atoms is lower resulting in a lower cp e - 
Fig. ^c) summarizes the discussion schematically and displays the link between the 
arriving particles, the erosion depth and the patterns formed. 
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Based on the data the three following conclusions are obtained, (i) Dot patterns 
form during 2keV Kr + ion beam erosion with *d = 30° for moderate co-sputtered Fe 
concentrations CF e and in the absence of reflected Kr particles, i.e. in the absence of 
a second component of energetic particles impinging from a different direction onto 
the substrate, (ii) Ripple patterns form on the Si substrate for high co-sputtered Fe 
concentrations: As the flux and the energy of the reflected Kr particles arriving in 
the ripple pattern area 1mm < x < 2.5 mm are small, their presence appears not to 
be necessary for ripple pattern formation. However Fig. [5] proves that the direction 
of impingement of the co-sputtered Fe atoms determines the orientation of the ripple 
wave vector k. Whether ripple orientation towards the stainless steel plate is due to the 
additional directed energy of the arriving Fe particles or due to the directed attachment 
of the Fe particles to elevations (ripple slopes) in their line of sight can not be resolved by 
the present experiment, (iii) Erosion by the direct beam together with strong sputtering 
from a different direction (here through the reflected energetic Kr particles) prevents 
formation of a pronounced and rough pattern. Although we detect a significant amount 
of Fe for x < 1 mm, the area remains rather smooth and the nanohole pattern is neither 
pronounced nor well ordered. 

As a global key result of our experiments we consider the fact that supply of a 
second chemical species is mandatory for pattern formation on Si in a large parameter 
range, specifically for d < 45°. Whether co-sputter deposition must be considered also 
for pattern formation of other materials or other energy ranges outside the 1 keV region 
as a hidden parameter remains to be investigated. To obtain a deeper understanding 
beyond this point and to disentangle the effects of a second chemical species and of 
additional energy from a second direction (at variance with the primary ion beam 
direction) dedicated experiments are necessary. It would be useful to design a situation 
where the second chemical species is deposited with thermal energies (e.g. through 
physical vapor deposition) simultaneously with a clean eroding ion beam. Such an 
experiment might also be able to establish whether the direction of the impinging co- 
deposited particles with respect to the ion beam and the sample surface is of relevance - 
irrespective of the energy of the co-deposited particles. As pattern formation is always 
the result of the interplay of a destabilizing mechanism and a stabilizing or smoothening 
one, the nature of the latter needs to be explored as well. We suggest to perform co- 
sputter deposition or co-evaporation experiments at very low (« 100 K) or elevated 
temperatures ~ 450 K to test the relevance of thermal diffusion. 

We believe that our experiments are relevant for future theoretical work. After 
the decade of continuum theory evolution there appears now to be a need for (i) a 
material parameter based description of pattern formation; (ii) for a theory that takes 
into account the absence of pattern formation for a large parameter space; (iii) a model 
which includes the effect of impurities for an explanation of pattern formation. 
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Figure 7. (a) Topography of the Si wafer measured by phase-shifting interferometry. 
Center of the long image axis identical with coordinate x. Dark blue contrast for the 
area with x < indicates the position of the steel plate and not height, (b) Large scale 
SEM image of the Si wafer. The position of the stainless steel plate during the erosion 
experiment is indicated, (c) Schematic sketch of the erosion situation with co sputter 
deposition (see text). 

5. Conclusions 

Room temperature ion beam erosion of Si through 2 keV Kr + shows no pattern formation 
for •& < 45°. For 60° < •& < 80° rough ion beam faceted ripple and roof tile patterns are 
formed. It is speculated that patterns observed for d < 45° are the result of unintentional 
co-sputter deposition of metallic impurities resulting from the use of broad beam ion 
sources with a divergent beam. The angular ranges of stability and instability as well 
as the observed faceting can be understood on the basis of the ^-dependence of the 
sputtering yield. 

Intentional co-sputter deposition of stainless steel from a plate hit together with 
the Si sample by the ion beam leads to a complex sequence of patterns in dependence of 
the distance from steel plate. These patterns form at angles where no patterns result 
for clean ion erosion. Nano hole, ripple and dot patterns are formed. Nano hole patterns 
are observed under the influence of strong additional sputtering by ions reflected from 
the steel plate and of co-sputtered steel. Ripple and dot patterns are formed in a 
distance range where reflected ions are largely absent. The ripple wave vector of the 
ripple patterns is found to be aligned to the direction of the impinging steel atoms. The 
formation of ripple and dot patterns is distinguished by the concentration of co-sputtered 
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material. Dots contain much less co-sputtered material than ripples. 

We believe that our results are relevant for future research in ion beam pattern 
formation on Si. They will foster additional experimental work to uncover the 
mechanisms of pattern formation due to the simultaneous deposition of a second 
chemical species. They also point to the need for well controlled erosion experiments 
and chemical analysis of the eroded surfaces to rule out impurity effects. Lastly, we 
hope they stimulate theoretical work, reconsidering the effect of the slope dependent 
sputtering yield for pattern formation and investigating the effects of impurities on 
pattern formation. 

Acknowledgments 

We thank Dietmar Hirsch for the SEM and SIMS measurements, Agnes Mill for optical 
micro phase interference measurements, Bernd Rauschenbach for support and useful 
discussions and Hubert Gnaser for his advice related to quantification of the SIMS data. 
We gratefully acknowledge the support by Deutsche Forschungsgemeinschaft through 
Forschergruppe 845. 

References 

[1] Wilson S J and Hutley M C 1982 Optica Acta 29 993 

[2] Flamm D, Frost F and Hirsch D 2001 Appl. Surf. Set. 179 95 

[3] Fouckhardt H, Steingoetter I, Brinkmann M, Hagemann M, Zarschizky H and Zschiedrich L 2007 

Advances in Opto Electronics Article ID 27316 1 
[4] Chaudhari P, Lacey J, Doyle J, Galligan E, Lien S C A, Callegari A, Hougham G, N D Lang 

PSA, John R, Yang K H, Lu M H, Cai C, Speidcll J, Purushothaman S, Ritsko J, Samant M, 

Stohr J, Nakagawa Y, Katoh Y, Saitoh Y, Sakai K, Satoh H, Odahara S, Nakano H, Nakagaki 

J and Shiota Y 2001 Nature 411 56 
[5] Zhang K, Rotter F, Uhrmacher M, Ronning C, Krauser J and Hofsass H 2007 New J. Phys. 9 29 
[6] Zhang K, Uhrmacher M, Hofsass H and Krauser J 2008 J. Appl. Phys. 103 083507 
[7] Liedke M O, Liedke B, Keller A, Hillebrands B, Miicklich A, Facsko S and Fassbender J 2007 Phys. 

Rev. B: Condens. Matter 75 220407 
[8] Teichert C, de Miguel J J and Bobek T 2009 J. Phys. Condens. Matter 21 224025 
[9] Kotai E 1994 Nucl. Instrum. Methods Phys. Res., Sect. B 85 588 
[10] Katardjiev I V, Carter G and Nobes M J 1989 Vacuum 39 1069 
[11] Carter G 2001 J. Phys. D: Appl. Phys. 34 Rl 

[12] Bradley R M and Harper JME 1988 J. Vac. Sci. Technol., A 6 2390 

[13] Madi C S, Davidovitch B, George H B, Norris S A, Brenner M P and Aziz M J 2008 Phys. Rev. 
Lett. 101 246102 

[14] Ziberi B, Cornejo M, Frost F and Rauschenbach B 2009 J. Phys. Condens. Matter 21 224003 
[15] Gago R, Vazquez L, Cucrno R, Varela M, Ballesteros C and Albella J M 2001 Appl. Phys. Lett. 
78 3316 

[16] Gago R, Vazquez L, Plantevin O, Sanchez-Garci'a J A, Varela M, Ballesteros M C, Albella J M 

and Metzger T H 2006 Phys. Rev. B: Condens. Matter 73 155414 
[17] Fan W B, Li W Q, Qi L J, Sun H T, Luo J, Zhao Y Y and Lu M 2005 Nanotechnology 16 1526 
[18] Fan W B, Ling L, Qi L J, Li W Q, Sun H T, Gu C X, Zhao Y Y and Lu M 2006 J. Phys. Condens. 

Matter 18 3367 



Is keV ion induced pattern formation on Si(OOl) caused by metal impurities? 15 

[19] Ziberi B, Frost F, Rauschenbach B and H6chc T 2005 Appl. Phys. Lett. 87 033113 
[20] Erlcbacher J, Aziz M J, Chason E, Sinclair M B and Floro J A 1999 Phys. Rev. Lett. 82 2330 
[21] Erlcbacher J, Aziz M J, Chason E, Sinclair M B and Floro J A 2000 Phys. Rev. Lett. 18 5800 
[22] Brown A D and Erlcbacher J 2005 Phys. Rev. B: Condens. Matter 72 075350 
[23] Gago R, Vazquez L, Cuerno R, Varela M, Ballesteros C and Albella J M 2002 Nanotechnology 13 
304 

[24] Habenicht S, Lieb K P, Koch J and Wieck A D 2002 Phys. Rev. B: Condens. Matter 65 115327 
[25] Keller A, Facsko S and Moller W 2008 New J. Phys. 10 063004 

[26] Keller A, Cuerno R, Facsko S and Moller W 2009 Phys. Rev. B: Condens. Matter 79 115437 
[27] Ziberi B, Frost F, H6chc T and Rauschenbach B 2005 Phys. Rev. B: Condens. Matter 72 235310 
[28] Ziberi B, Frost F, Hoche T and Rauschenbach B 2006 Vacuum 81 155 
[29] Ziberi B, Frost F and Rauschenbach B 2006 J. Vac. Sci. Technol., A 24 1344 
[30] Makeev M A, Cuerno R and Barabasi A L 2002 Nucl. Instrum. Methods Phys. Res., Sect. B 197 
185 

[31] Facsko S, Bobck T, Stahl A and Kurz H 2004 Phys. Rev. B: Condens. Matter 69 153412 

[32] Castro M, Cuerno R, Vazquez L and Gago R 2005 Phys. Rev. Lett. 94 016102 

[33] Vogel S and Linz S J 2007 Phys. Rev. B: Condens. Matter 75 085425 

[34] Chan W L and Chason E 2007 J. Appl. Phys. 101 121301 

[35] Wehner G K and Hajicek D J 1971 J. Appl. Phys. 42 1145 

[36] Punzel J and Hauffc W 1972 Phys. Status Solidi A 14 K97 

[37] Robinson R S and Rossnagel S M 1982 J. Vac. Sci. Technol. 21 790 

[38] Rossnagel S M and Robinson R S 1982 J. Vac. Sci. Technol. 20 195 

[39] Tanemura M, Yamauchi H, Yamanc Y, Ohita T and Tanemura S 2004 Nucl. Instrum. Methods 

Phys. Res., Sect. B 215 137 
[40] Zhang K, Hofsass H, Rotter F, Uhrmacher M, Ronning C and Krauser J 2009 Surf. Coat. Technol. 

203 2395 

[41] Ozaydin G, Oczan A S, Wang Y, Ludwig K F, Zhou H, Hcadrick R L and Siddons D P 2005 Appl. 

Phys. Lett. 87 163104 
[42] Ozaydin-Incc G and Ludwig Jr K F 2009 J. Phys. Condens. Matter 21 224008 
[43] Teichcrt C, Hofer C and Hlawacek G 2006 Adv. Eng. Mater. 8 1057 

[44] Sanchez-Garcfa J A, Vazquez L, Gago R, Redondo-Cubero A, Albella J M and Czigany Z 2008 

Nanotechnology 19 224009 
[45] Sanchez-Garcfa J A, Gago R, Caillard R, Redondo-Cubero A, Martin-Gago J A, Palomares F J, 

Fernandez M and Vazquez L 2009 J. Phys. Condens. Matter 21 224009 
[46] Polop C, Rosiepen C, Bleikamp S, Drese R, Mayer J, Dimyati A and Michely T 2007 New J. Phys. 

9 1 

[47] Horcas I, Fernandez R, Gomez-Rodriguez J M, Colchero J, Gomez-Herrero J and Baro A M 2007 

Rev. Sci. Instrum. 78 013705 
[48] Eckstein W 1994 Radiat. Eff. Defects Solids 130 239 
[49] Ducommun J P, Cantagrel M and Moulin M 1975 J. Mater. Sci. 10 52 
[50] Carter G, Colligon J S and Nobes M J 1973 J. Mat. Sci. 8 1473 
[51] Carter G, Colligon J S and Nobes M J 1977 Radiat. Eff. Defects Solids 31 65 
[52] Wei Q, Lian J, Boatncr L A, Wang L M and Ewing R C 2009 Phys. Rev. B: Condens. Matter 80 

085413 

[53] Carter G, Nobes M J and Katardjiev I V 1992 Philos. Mag. B 66 419 
[54] Stevic F A and Wilson R G 1991 J. Vac. Sci. Technol, A 9 3064 
[55] Wilson R G 1995 Int. J. Mass Spectrom. Ion Processes 143 43 



